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Abstract 

We present a flat {K = 0) cosmological model, described by a perfect fluid with the "con- 
stants" G, c and A varying with cosmological time t. We introduce Planck's "constant" h 
in the field equations through the equation of state for the energy density of radiation. We 
I then determine the behaviour of the "constants" by using the zero divergence of the second 



member of the modified Einstein's field equations i.e. div{^T- + d^A) = 0, together with the 



PACS numbers: 



. equation of state and the Einstein cosmological equations. Assuming realistic physical and 

0^ \ mathematical conditions we obtain a consistent result with he = constant. In this way we ob- 

tain gauge invariance for the Schrodinger equation and the behaviour of the remaining constants" . 
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I. INTRODUCTION. 



We study the behavior of the G, c and A constants in the frame described by a model with 
an energy-momentum tensor defined by a perfect fluid, considering the symmetry implied 
by a Friedman-Robertson- Walker metric (we restrict to the flat 3-space case, = 0, as has 
been recently observed in the CMBR see' ). 

In order to incorporate in the field equations constants like Planck's h, we consider for 
the radiation-dominated case an energy density given by the equation of state 

p = aO^ (1) 

where a represents the radiation constant and 6 the temperature. Using this expression for 
the energy density we can get a possible time variation for h from the field equations. 

We use the modified Einstein's field equations (see-^, ' and^) considering the possibility 
that all constants may vary with time t, and only with time as inferred from the homogeneity 
and isotropy of the metric. We then impose a zero value for the covariant divergence of the 
second member of the field equations 

div {Ti + 5iK) = (2) 

This equation, together with the equation of state and the field equations, will give us the 
behavior of the "constants" . To solve these equations we impose some restrictions and this 
will drive us to various possibilities. And we do not consider starting with the condition 
of zero covariant divergence for the energy- momentum tensor. This condition has been 
analyzed (in a phenomeno logical way) by various authors (see*^) that explain in this way the 
creation of matter and entropy as due to the time-variation of G, c and A only. 

Our purpose is not to study thermodynamic properties of the model, but to determine 
it including h in the field equations. We impose two physical and mathematical realistic 
conditions: 

• It will always be true the G and c relation 

^ = const. (3) 

which is a general covariance principle. We know that this condition is verified (not 
imposed) in the case of a bulk viscous fluid. And we restrict ourselves to the only case 
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such that the parameter controlhng the viscosity describes a behavior topologically 
equivalent to the case given by the classical FRW applied to a perfect fluid (our case) 
(see^ and"'). 

We impose the following relation for the cosmological "constant" obtained through 
the similarity group that applies to the equations (a mathematical condition, as a 
self-similarity relation) 

(4) 

Another option that we have, which is equivalent to the above as we will see, and 
already analyzed in the literature, is the following one (see'): 

A oc (5) 



The paper is organized as follows: In section 2 we present the equations that define 
the model and describe the different simplifying hypothesis to be imposed in the subsequent 
sections. In section 3 we study a model without the condition of zero covariant divergence for 
the energy-momentum tensor. We will see that with the results obtained the fine structure 
constant a continues being constant (see'^), despite the time- variation of all the constants 
in it. In section 4 we will also see that Schrodinger and Maxwell's equations remain gauge 
invariant (see") with the time varying functions found for the constants in our study. Also 
our model does not have the Planck's problem. 

II. THE MODEL. 
A. The equations: 

We will use the field equations in the form: 

Rij - \9^jR = ^^^^^Tij + A{t)gij (6) 
Where the energy momentum tensor is: 

Tij = {p + p) UiUj - pQij (7) 
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where p = up in such a way that u G [0, 1] . 
The cosmological equations are now: 

where / stands for the scale factor. Applying the covariance divergence to the second 
member of equation (6) we get: 

=(^- ^] - ^^^^^^^^ (W) 

'■'^ \ c G J ' SnG ^ ^ 

that simplifies to: 

p' + 3iuj + l)pH = - ^ - + 4p- (11) 

T^. ' 

^■,J 

f 

where H stands for the Hubble parameter {H = —). The last equation may be written in 
the form: 

p' + 3(c. + l)pi/+^ + p^-4p- = (12) 

QTXLx Lx C 

or the equivalent 

^ + S{u + 1)H+^ + ^-A'- = (13) 
p SvrGp G c 

In this equation we can take a general energy density p or use the black body equation 
of state p = a9^ where a = 5—^ so that equation (13) is now: 



9' 

4- -3 



c h 



'^1 J 



+ 3(^ + 1)^+^-777747^7 + 77-4- = (14) 



M^Gk%e^ G 



c 



B. Simplifying hypothesis 

1. Hypothesis number one (HI). 
We introduce 



G 

— = const. = B (15) 



c 
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with [B] = LM-^. In this way we obtain directly a term that controls the behavior of the 
energy density. 

P=^ 6 G M (16) 
If we take p = — ^ and introduce p = aO'^ we arrive at 



b 



1/4 



a9^^e= i-] t-/^ (17) 



V Ba 

2. Hypothesis number two (H2) 

We can make a previous hypothesis on the behavior of the Planck's constant as follows: 

• ch = const, or equivalently h = — , where A is a proportionality constant. We are 

c 

motivated to this relation because it leaves the radiation constant a as constant with 
a constant ks = const, (see') 

• h = Fc where [F] = L^M^. This hypothesis is introduced by some author to explain 
a possible variation of the fine structure constant. We will see that such hypothesis is 
irreconciliable with our model fsee"). 



3. Hypothesis number three (H3) 

As a mathematical possibility (self-similar relation), obtained by means of the similarity 
group acting on the equations, we will assume the following behavior for the cosmological 
"constant" : 

'''' 

where ci G M. This hypothesis is a strong one since in some way it forces to the relation 
f{t) = c{t)t, but this is a desired relation to get the causality principle preserved. 
With this hypothesis we now deal with the equations 

III. MODEL WITH div{T) / 0. 



As stated in the introduction this type of model has been recently studied by different 
authors. The general idea followed by Harko et al. (see'') is to study the consequences of 
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the time variation of the constants in order to solve some problems of the Standard Hot 
Big-Bang Model (SHBB). 
By imposing the condition 

dw{T^) > (19) 

already used by Rastall (see^°) as of today there has been no proof that the equality div(T) = 
is verified experimentally, hence they impose the condition div{T-) > a where a is a 
certain positive function. Rastall imposes a law div{T-) = XRi i.e. proportional to the 
gradient of the scalar curvature. Hark et al. interpret the time variation of the constants 
as responsible for matter creation, and therefore for entropy, working with the following 
differential inequality 



^i = -^-p- + 4p-:>0 (20) 



A'c^ G' , c' 

But in order to solve it they have to impose a previous behavior for the "constants" they 
analyze. 

We will follow a different approach here (see^'^). We will impose that the covariant 
divergence of the second member of the Einstein's equations be zero, i.e. 

^ + 3(c. + l)i7+-^ + ^-4- = (21) 
p SnGp G c 

and in order to solve it we impose the hypothesis already cited. 



G 

A. Hypothesis one = const. = B together with hypothesis two A 



C2 ^ -""^ C^t)f^ 

In this case the energy density is of a general nature and the equation to be solved is 
then (13) 

^ + 3(^ + l)if+^ + ^-4^ = (22) 
p SnGp G c 

which we now analyze. 

Smce — = const. = B then p = — - so that — = — 
c2 ^ Bt^ P t 



Taking into account that ^ = 2— we get the following expression 



A'c^ d 



SnGp Anb 

Using the field equation 

f2 ^2 



c' i 

-c^-t 



(23) 



3^^ = ^p + c'A (24) 
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we obtain H and therefore /. 



3H' 



Hence we get / = K^t^ with x 
equation (22) it is obtained: 

-2 



MBh d 
1 

87r6 + d\2 
d c' 



t 



(25) 



And rearranging all these results in the 



, X a c d \ d 
- + 3 + 1 - - — - - — - + 2- 
t t Anb c Anh t c 



(26) 



1 


' d 


d 




^ + 2 




1 ~ 


A-Kh 


c 



simplifying: 



c 
c 



l2Tih{uj + l)x - %TTh - d 



Sub + d 
c = K^t^ 



(27) 



(28) 



(29) 



integrating it we obtain easily: 

^ \127rb(uj + l)x-87Tb-d' 

where ^ = g , , , 

L Sirb + d 

Another way to look at it is as follows: since we are working with the constant B when 
getting / we have obtained another constant i^^ with dimensions [i^^] = LT~^. With these 
two dimensional constants we have the solution without the need of integrating. We can 
check it with the following result obtained by means of gauge relations. 

The group of governing quantities is 971 = OJt {B, K^, t}, and the results are: 

c oc KJ''^^ 
G oc 

p oc B~H-'^ 
f = KJ- 

A oc 



(30) 



so that we can check the following relations f = d:, G = Bc^. We also see that A oc /" 

c"^t~^. With the result x = - we obtain: 

2 



c oc 



G oc t-^ 



f oc t^/^ p oc t" 



A oc 



(31) 
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Finally we see that we can regain the law div{T) = taking into account that with / oc t^^^ 
and p oc the relation pf^ = const, is satisfied, a relation obtained integrating the equation 
div{T) = 0. Our results are compatible with this relation. 



B. Hypothesis 3 A 



together with — = const. = B and p = aO^ 



c^{t)t^ 



Since the constant B has dimensions [B] = LM ^ we can get the dimensionless monomia 
pBt^ 

TTi = — - — where 6 G M. With these hypothesis and tti equation ((14)) simplifies to: 



9' 

4- -3 



15dc^[c't + c]h^ G' d 



(32) 



that has no immediate integration. We have to take into account the field equations (9). 



3if^ = 87r6t-2 + dr 



(33) 



from this we get / = KJ:'^ where x = ( — j ^ and substituting in (32) together with 



G = Bc^ we get 



4- 



3{uj + l)x _ ISdc^ [c't + c]h^ _ d _ 
^ t 47r3 Bk%eH^ 7~ 



(34) 



i.e. one equation with 3 unknowns. Similarly to the previous case we can take into account 
the following group of governing quantities 071 = 971 {S, i^^, t} in such a way that 



c oc KJ'^-^ 



(35) 



etc (same solution as the previous section, see 30) 

If we want to integrate equation (34) we have to take a decision on the behavior of he 
constant h. 



A h! d 
1. Taking dh, = const. ^ h = — then — = yielding: 

c he 

9' 3{u + l)x 15dA^dt + c] 
^1 ^ i 47r3 cBk%dH^ 



2^ = 

c 



(36) 



Also if p = aO^ and p 



Bt^ 



Bt^ 



aO^ we have: 



TT^Bt^ 



■B J 



(37) 



And substituting into the previous equation we get : 

-2 3(cu + l)x 15d [c't + c] 



and simphfying 



t t 



-2 3cu + lx 

— + ^ — 

t t 



47r3 ^ A5A% 



(38) 



d 

Airb 







(39) 



we see that this is the same equations as (27 ) and therefore will have the same solution 



We can consider another possibility. Take the group of governing quantities 
^)Jl = DJl {K^, A, t} where is the proportionality constant obtained from / = K^t"^ 
and A is the constant establishing the relation between h and c. The results obtained 
by means of the gauge relations are: 

G oc K^A-H^'^-'^ 
c oc K^t''-^ 
h oc K~^At^~'' 

kse oc K-^M->^ ^^^^ 
p oc K-'^At"^'^ 

A oc K~H~^^ 
e^e^^ oc A 

9 
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Where mj comes from the energy density definition pE = — -p; — [n stands for the 

■' 

particles number) and c^Eq from the definition of the fine structure constant a. We 

G 1 

can check that we recover the general covariance property — = const, iif x = -. 

c 2 

Similarly we can see that the following relations are satisfied: p = aO^^ p = Aj~^ 

(equivalent to div{T) = 0), A oc and f = ct (no horizon problem). And finally 

e^Sn^ oc const. In this way the fine structure constant a oc = const., is a true 

eoch 

constant (see^). (For a more detailed analysis on the behavior of the electromagnetic 
constants see' and"). With the value x = - we get 

cocri/2, hoct^/\ Goct-\ kB0(xt-^/^ , ^ 

(41) 

fo^t^/'^, poct^^ mi oct^/2 
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i.e. same result as in (31). We emphasize that all these results coincide with the work 
of Midy et al (see'^). 

Another alternative is to consider the governing quantities DJl = DJl {B, A, t} such that 

[A] = [A^] with u = -. So that we get the same problem as with the condition 

3 

div{T) = ( Note p = A^f-^^^+^^) This would imply no previous hypothesis on the 
behavior of h. 



2. Assuming h = Fc the equations (34) is now: 



9' 

4--3 



3(cj + l)x 15d \c't + c] c' 
-I !^ i 1 i 2— = 

t 47r3 Bk%eH^ c 



(42) 



With the relation p = aO'^ and p 



Bt^ 



a9^ 



1 1 



So that 



we obtain: 



9' c' 2 

4- = 6 

e c t 



'■Bt^ 



-2 3(c<; + l)x d 



2- = 

c 



(43) 



(44) 



(45) 



t t A-Kh 

Same solution as before (39). We are getting the same equation, and the same solution, 
because we are imposing p oc t~^. 

If we check this hypothesis with the set of governing quantities OJl = 071{-B,-F, t} we 
get: 

G oc B'^Ft-'^ 
c oc B^I^F^IH-^ 
f = const, 
p oc B-H-^ 
ksO oc B^/^F^/H~^ 
h oc B^'^F'^'H-^ 



(46) 



As we see these results are discouraging. In the previous results all the hypothesis were 
compatible. Here they are not. 
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IV. GAUGE INVARIANCE. 



We will check the gauge invariance of the Schrodinger, Klein-Gordon and Maxwell equa- 
tions fsee^^). 



1. Consider Schrodinger equation: 



2m 2n 



(47) 



2. Klein-Gordon equation: 



1 d'^ip rn^c^ 



+ -^4^ = 



3. Maxwell equations are ( see^'^): 



(48) 



rotB = —dtE rotE = -dtB 



(49) 



divB = divE = ^ 

So 



(50) 



(One of us (A,A-F) have deduced a modified Maxwell equations that remains gauge 

invariant too, for a more detailed study see^^) 

These equations are all independent from one to another. By introducing a characteristic 
length L and time T we can write all these equation in a dimensionless way. 



Schrodinger equation is then : 

/ = LC 



V = -(5 t = TT 

1j 



(51) 



and therefore: 



Klein- Gordon equation 



JJ 



h"^ h 



2mL2 



5^ij 1 ^ 



(52) 
(53) 

(54) 
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While Maxwell's equations give 

B = B*P E = E*e 



Q 



(55) 



R* TP* f)(: 

The invariance of all these equations is ensured if: 

hT T h 

1. — — = const, i.e. ~ — 

h r, h 

2. L « cT, mc « — , mc ^ — 

E* Q 

3. L cT and — ~ A relation that is consistent with the definition of an 

L EqR-^ 

electric field due to an electric charge Q. 
As we see these are the results obtained elsewhere (41) i.e. 

/ oc m^h(x (58) 

We can check now that our model does not have the so called Planck's problem. The 
system behaves now as: 

And since the radius of the Universe / at Planck's epoch coincides with the Planck's length 

f{tp) ~ Ip (60) 
and the energy density at Planck's epoch coincides with Planck's energy density 

p{t,) ^p,^ (61) 

where Pp = mpC^/lp, we have no Planck's problem here. 
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V. CONCLUSIONS. 



We have presented here a flat cosmological model in which many "constants" are consid- 
ered time-varying. We have found that the only possibility compatible with the established 
hypothesis is that the product ch be a constant. In this way the equations of physics are 
gauge invariant. The case considered here is one of a radiation-dominated epoch. However 
the time variations found are considered to be of some general value in the sense that an 
epoch that is not radiation dominated may conserve some of the properties found. Further 
analysis proves that this is the case and that the result / = t^^^ has a general validity. 

Despite of starting with conditions different from div{T) = we have found that this 
relation is verified for the quantities of our model. With these results we can say that we 
can consider adiabatic creation mechanisms of matter in order to get rid of the entropy 
problem. In our model the main quantities present a similar behavior as the FRW classical 
models, but we have the advantage of no such problems as the horizon one. The behavior 
found for the fine structure constant is that it is a true constant, and that we can get rid 
of the Planck's problem. The horizon problem is solved as soon as the relation / = ct is 
satisfied. Therefore the Universe expands at the speed of light. 
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